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□ Fundamentals Of GIS □ 



OVERVIEW OF GIS 



G I S 

Geographic Information Systems 

An INTEGRATED SYSTEM of COMPUTER HARDWARE and SOFTWARE 
coupled with PROCEDURES and a HUMAN ANALYST which together 
support the CAPTURE, MANAGEMENT, MANIPULATION, ANALYSIS . 
MODELLING, and DISPLAY of SPATIALLY REFERENCED DATA 



GEOGRAPHIC INFORMATION SYSTEM 



DATA 
ENTRY 



DATA 
STORAGE 



DATA 
MANIPULATION 




MAPS DIGITAL TABULAR 
DATA SUMMARIES 



'ANSWERS" 
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G I S Capabilities: 

□ Query for Location 

"Show me all the countries of South America that have a 
population greater than 20,000,000. " 

□ Query for Condition 

"Display the population of each country I point to on the 
map." 

□ Trend Analysis 

"Show me where the census blocks are that have 
experienced more than a 50% population change 
between 1980 and 1990." 



□ Patterns Analysis 

"Calculate the fragmentation index for all the forest 
patches in the municipio." 

□ Modelling 

"Which route for the new highway has the lowest cost in 
terms of losses of housing, prime farmland, and wetlands, 
while minimizing the needs for cutting and filling." 
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Query: Population > 20,000,000 



Simple Query 



l£\ South America 



menc 



Coitfitw Name 



^JnJilJ 



Polygon 


j Argentina 


j 33796870! 


A 


Polygon 


Brazil 


! 151525400 




Polygon 


j Colombia 


[ 34414590 




Polygon 


j Peru 


j 24496400 




Polygon 


j Bolivia 


7648315 




Polygon 


| Chile 


j 13772710! 
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d Databas 



action 



l£\ South America 



^ Attributes of South America 


1- n 


-I 




F'opiriaf/as? 




Polygon j Argentina 


33796870! 




Polygon j Bolivia 


7648315 




Polygon i Brazil 


151525400! 




Polygon j Chile 


137727101 




Polygon j Colombia 


34414580! 
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Point 



Well or Town 



Line 




Stream or Road 



Area 




City or Field 



Volume 





>4 



















Fertilizer or Yield 



Geometric Classes of Data 



DATA TYPES 



VECTOR 



RASTER 



Point = Position, no area 



Point = 1 cell 



Line = Length, no width 



Line = Multiple cells joined at 

edges or corners, usually 
with only 1 or 2 neighbors 



Polygon = Area and perimeter 



Polygon = Group of contiguous 

cells joined at edges or 
corners 
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VECTOR STRUCTURE 
Advantages 

□ Good representation of the landscape being mapped 

□ Topology can be completely described, including network 
linkages 

□ Great looking graphics ("Looks like a map is supposed to") 

□ Generalization of the graphics is possible while still maintaining 
the great look ("What the map reader doesn't know won't hurt 
them") 



RASTER STRUCTURE 
Advantages 

□ Overlaying maps is easy and "perfect" (i.e. no possibility of 
sliver polygons developing since all raster cell borders are 
coincident 

□ Integration of remotely sensed imagery (satellite images or 
scanned airphotos) is straight-forward 

□ A huge variety of complex spatial analyses are supported 

□ Software is generally cheaper and easier to learn compared to 
vector GISs 
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TOPOLOGY 

Geometrical relationships between spatial objects (Points, Lines, and 
Areas), such as adjacency, that are not altered by distortion, as long as the 

surface is not torn 

Example of "Built" Topology (from Arc/Info) 




Node D 

Arc ID 

Polygon ID 

arc diqit'zea 
in direction 
of arrow 



— w — 


Arc 
ID 


Left 
Poly 


Right 
Poly 


From 
Node 


To 

Node 








1 


A 


0 


c 


a 




2 


A 


0 


b 


c 




3 


C 


A 


b 


a 




4 


0 


C 


d 


a 




5 


C 


B 


d 


b 
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B 


D 
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e 


Polygon 
ID 


Number 
of arcs 


List 
of arcs 


B 


0 


d 


c 


A 


3 


-1,-2,3 




B 


4 


2,-7,5,0,-6 


C 


3 


-3,-5,4 


D 


1 
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Vector GIS ATTRIBUTE DATABASE 



ID 


Name 


Quality 




0 


N/A 






1 


Fir Lake 


High 





ID 



ID 


Species 


Size 




0 


N/A 






1 


Jack Pine 


Pole 




2 


Spruce 


Sapling 




3 


Jack Pine 


Sapling 





Type 



N/A 



Depth 



Lakes 




Forest Types 





1 t 


1^ 










2 











Soils 




Type 


Drainage 




A 


Fair 




B 


Poor 




C 


Good 




R 


Rock 





Related Table 
"Type" is 
KEY FIELD 
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ENTITY RELATIONSHIPS 

Point - Point Is nearest to ... 

Interacts with ... 

Point - Line Is on ... 

Is nearest to ... 

Point - Area Is within ... 

Is adjacent to ... 

Line - Line Intersects ... 

Is upstream of ... 

Line - Area Crosses ... 

Is contained within 
Is nearest to ... 

Area - Area Is adjacent to ... 

Overlaps ... 

Is enclosed by ... 
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A Classification of GIS Functions 

□ Analysis of Spatial and Attribute Data 

- Non-spatial analyses 

Attribute query and display 
Map retrieval and display 
Attribute classification 

Map measurements (distance, direction, area, etc.) 

- Spatial analyses 

Overlay operations 
Neighborhood functions 

Distance and Connectivity functions 
Contiguity measures 
Proximity analysis 
Network analysis 
Spread functions 
Seek operands 
Intervisibility analysis 
Solar illumination calculation 
Perspective view 

Search operations 

Line-in-polygon; Point-in-polygon 
Topographic functions 
Thiessen polygons 
Interpolation 
Contour generation 



This classification has been adapted from: 

Aronoff, Stan. 1989. Geographic Information Systems: A Management 
Approach. Ottawa, Ontario, Canada: WDL Publications. 294p. 
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A Classification of GIS Functions 



□ Maintenance of Attribute Data 

Format conversions 
Database error checking 
Database editing 



□ Maintenance of Spatial Data 

Format conversions 

Geometric transformations 

Projection conversions 

Conflation 

Edge matching 

Editing of graphic elements 

Line coordinate thinning 



□ Output functions 

Map annotation 
Text labels 

Texture patterns and line styles 

Graphic symbols 

Plotting 

Printing (laser printers, color inkjet printers, etc.) 



This classification has been adapted from: 

Aronoff, Stan. 1989. Geographic Information Systems: A Management 
Approach. Ottawa, Ontario, Canada: WDL Publications. 294p. 
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Spatial Analyses 



Basic Functional Classes 



Reclassifying Maps 



j 



Vector and Raster 



□ Overlaying Maps 



□ Measuring Distance and Connectivity 



□ Characterizing Neighborhoods 



Position 



Vector and Raster 



only those in the NW" 




Value 



Vector and Raster 



change feet to meters" 



elevations bet-ween 20 & 40 feet" 
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RECLASSIFYING MAPS 



Size 



Vector and Raster 



larger than 80 acres' 



^7/ 




Contiguity 



Raster Only 



"work with individual 
members raihers than 
the class as a whole" 




Vector and Raster 



Shape 

• Spatial Integrity 




[# Holes - (# Fragments - 1 )] 



Boundary 
Configuration 




smallest ^/-^^ large 
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SPATIAL ANALYSES 



Basic Functional Classes 



□ Reclassifying Maps 



Overlaying Maps 



j 



Vector and Raster 



□ Measuring Distance and Connectivity 



□ Characterizing Neighborhoods 



Point by Point 



Vector and Raster 



n 



93 



piercing — needle 

approach ^ 



All locations in the coverage or grid 
are evaluated. The results extend 
to the spatial limits of the input maps. 
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OVERLAYING 



Region-wide 



Vector and Raster 



7> 



cookie— cutter 
approach 



All locations in the 
coverage or grid 

are not necessarily 
evaluated. The results 
are constrained to the 

spatial nature of the 
v reference map. , 





Topologic 
Overlay 



"co — oc currence 



Vector Only 



jnappvng 




No"te "the new NDDE - "topaloglcalty 
correct areas were created 
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V Forests 



_ n x 





R hre r 








L 


Forest Stands 






| | Harvestable 






| | Protected 











Intersect - Query: Forests = Harvestable and Erosion = Low 



n x 



; rv * r 



^ RK, 

A/ 

i/f Inte rse ct.s hp 



Forest Stands 
| | Harvestable 
I I Protected 



•2 Intersect 













Polygon 


7 


High 


Protected 




Polygon 


2 


Low 


Harvestable 




Polygon 


2 


Low 


Harvestable 




Polygon 


8 


Low 


Protected 




Polygon 


9 


Low 


Protected 


_ 


Polyoon 


10 


Low 


Harvestable 




< ITT 
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OVERLAYING MAPS 



Boolean 
Overlay 



Vector and Raster 



Intersection 

"Forest AND Soil A" 




Union 

"Forest OR Soil A' 



yyyyyyyyyyy 

^££y^yy 



•Negate 

"Forest, but NOT on 
Soil A" 

•Exclusive Or 

"Forest XOR Soil A" 




Note: Most relational database management systems support a "SWAP" 
function which selects the currently unselected items in the database. In 
the example above, the "SWAP" function would return the white areas after 
the hatched areas had initially been selected. 
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SPATIAL ANALYSES 



Basic Functional Classes 



□ Reclassifying Maps 



□ Overlaying Maps 



□ Measuring Distance and Connectivity 



) 



□ Characterizing Neighborhoods 



DISTANCE MEASURES 



Point to Point C 



'ector and Raster 



"How far is it from A to B?" 

USEFUL. BUT LIMITED 



Proximity 

How far from the forest is every 
location on the map?" 

POWERFUL AND VERY DIFFICULT TO DO BY HAND 



Raster Only 
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DISTANCE MEASURES 



Movement 



Raster Only 





Start Location Travel Time 

(Impedance = 0) 




Movement 



Raster Only 




A = gently sloping grassland 
B = steeply sloping forest 
C " gently sloping forest 



1 = easy hiking 

2 = harder hiking 

3 = difficult hiking 




Start Location 



Travel Time 

(with variable impedances) 
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Least-Cost Pathway Analyses 

□ Create a FRICTION map from one or more 
existing coverages 

a Create a COST SURFACE map by executing 
the PROXIMITY analysis WEIGHTED BY the 
FRICTION map 

□ Execute the LEAST-COST PATHWAY module 
from one or more starter entities (points, 
lines or areas) to a destination. It will FIND 
the one (or more) routes that ACCUMULATE 
the least cost. 



c) 1999 David P. Lusch, PhD O Center For Remote Sensing & GIS O Michigan State University 23 



□ Fundamentals Of GIS □ 




Base cost 



Base cost x 4 



Base cost x 5 



Base costx 1000 



Friction Surface 

(reclassed from the 
Land Use map) 
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Cost Surface 
(Proximity x Friction) 



* 
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Least-Cost Pathway 




Least-Cost Pathway 
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CONNECTIVITY 



Buffering 



Vector and Raster 





establish limited use zones 



100 m 

200 771 

500 m 



1 st— order streams 
2nd— order streams 
3rd— order streams 



Q Buffers 



^ Parana ^ 

A/ I 

^ 100 Miles 



1MB 



^ Buffers 



VT Parana 

A/ 

^ Miles 

| | 100 
i | 200 
300 
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Point j 


1 


Iron Ore 


0.645 


Parana 


888621 


Point j 


2 


Gravel 


1.486 


Parana 


888621 


Point j 


3 


Gravel 


2.460 


Parana 


888621 


Point j 


4 


Iron Ore 


1.821 


Parana 


888621 


Point j 


5 


Gravel 


1.615 


Parana 


888621 


Point j 


6 


Gravel 


1.835 


Parana 


888621 




<l 
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SPATIAL ANALYSES 



Basic Functional Classes 



□ Reclassifying Maps 



□ Overlaying Maps 



□ Measuring Distance and Connectivity 



Characterizing Neighborhoods 



Raster Only 



CHARACTERIZING 
NEIGHBORHOODS 




Neighborhood = 
User— Defined Area 
of Scrutiny 

"Roving Window" 
(e.g. ... per square mile; ... per hectare) 
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WINDOW OPERATIONS 



Raster Only 



Slope 

Slope Aspect 

Maximum, Minimum 

Mean, Median, Mode 

Standard Deviation 

Majority, Minority 

Total Count 

Diversity 

Large number of Spatial Pattern 
or Texture Indexes 
(e.g. Dominance Index, Relative Richness, etc.) 
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s© x Neighborhood Functions 



nJU 



h|Ti 



1489 - 
14192 - 
1495 - 
1498 - 
1501 - 
1504- 
1506 - 
1509 - 
1512 - 
1515 - 
1518 - 
1521 - 
1524- 
1526 - 
1529 - 



1492 
1495 
1498 
1501 
1504 
1506 
1509 
1512 
1515 
1518 
1521 
1524 
1526 
1529 
1532 




^ Neighborhood Functions 



si 



i/f As p e ct of D e m 
| n Flat (-1) 



North 

N orth east 

East 

So uth east 

South 

So uth west 

West 

N orth west 



^ Slope of Dem 

| H 0.1 - 3.2 
| | 3.2 - 6.2 
| | 6.2 - 9.3 

| j 9.3 - 12.4 
| | 12.4 - 15.4 
| | 15.4 - 18.5 
I I 18.5 - 21.6 



Aspect 



Slope 
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Examples of Neighborhood Operations to Determine 

Spatial Pattern 

Each of the following measures are calculated within a 3 x 3 pixel window which systematically roves 
throught the data set. The outcome calculation is assigned to the center cell in the window (in the output 
file), then the window moves over one pixel along a row and recalculates a new value for that center- 
pixel location, etc. 



Relative Richness = n / (nmax) x 100 

where n = number of different classes present 

Diversity = -sum[(p) x ln(p)] 

where sum = the sum over all classes; p = proportion of the footprint in each class; In = natural logarithm 

Dominance = Hmax - H 

where H = Diversity; Hmax = maximum diversity = ln(n); n = number of different classes present; 
In = natural logarithm 

from: Turner, M.G. 1989. Landscape Ecology: The Effect of Pattern on Process, Anttu, Rev. Ecol. Syst, 20, 171-197. 
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